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Abstract. Allowing the energy of a gravitational field to 
serve partially as its own source allows gravitating bod- 
ies to exhibit stronger fields, as if they were more mas- 
sive. Depending on degree of compaction of the body, the 
field could be one to five times larger than the newto- 
nian field. This is a comfortable range of increase in field 
strength and may prove to be of convenience in the study 
of velocity curves of spirals, of velocity dispersions in clus- 
ters of galaxies and in interpreting the Tully-Fisher or 
Faber-Jackson relations in galaxies or systems of galax- 
ies. The revised gravitation admits of superposition princi- 
ple but only approximately in systems whose components 
are widely separated. The revised dynamics admits of the 
equivalence principle in that, the effective force acting on 
a test particle is derived from a potential, and could be 
eliminated in a freely falling frame of reference. 
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The only justification for our concepts and system of 
concepts is that they serve to represent the complex of 
our experiences. Beyond this they have no legitimacy. I 
am convinced that the philosophers have had a harmful 
effect upon the progress of scientific thinking in remov- 
ing certain fundamental concepts from the domain of 
empiricism, where they are under our control, to the in- 
tangible heights of the a priori. Albert Einstein, 1922. 

1. Introduction 

The problem of dark matter surrounding spiral galaxies, 
wrote Van Albada et al. (1985), is one of the most enig- 
matic questions in present day astrophysics. A number of 
years of intensive research have brought little or no clari- 
fication. 

Earlier suggestions to resolve mass discrepancy were 
on the conservative side. Bahcall and Cassertano (1985) 
spoke of missing mass as synonymous to missing light. 
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Tremaine and Lee (1987)with dark matter meant the mat- 
ter whose existence was inferred only through its gravita- 
tional effects. Critics were vocal. Outstanding among them 
Milgrom (1983 a,b,c), the architect of MOND, spoke of 
the dark side of the dark matter hypothesis and of its ar- 
bitrariness in that one invokes the dark matter in the cor- 
rect amount and spatial distribution needed to explain the 
mass discrepancy in each and every case for itself (1987). 

In the course of the past two decades acurate ob- 
servations of velocity curves of spiral galaxies have be- 
come available. See, e. g., Begeman et al. (1991), Sanders 
(1996), Sanders and Verheijen (1998), Mc Gauph and de 
Blok (1998). Credible velocity dispersions of galactic sys- 
tems have emerged. See, for example, Tullly et al. (1996) 
for data on Ursa Majoris cluster. On the theoretical side 
many ingenious and bold candidates for dark matter have 
been offered. A list compiled by Ostriker and Stcinhardt 
(2003) includes items such as cold collisionless dark mat- 
ter, strongly self interacting dark matter, warm dark mat- 
ter, repulsive dark matter, fuzzy dark matter, self annihi- 
lating dark matter, decaying dark matter, etc. 

In spite of all these developments, however, there is 
no clue as how to detect the dark matter and get to the 
physics that it obeys. The issue does not seem less enig- 
matic than what von Albada and his colleagues described 
in 1985. Arguments pro and con are the same as a quar- 
ter of century ago. Proponents of dark matter assume the 
validity of Newtonian dynamics at all distance scales and 
look for one or another form of hypothetical matter to 
provide the missing gravity. Skeptics, on the other hand, 
see no logic in resorting to a concept that solves the riddle 
of dynamics but remains unamenable to further validation 
by any other known physical means. 

A historical reminder might be timely. In the closing 
decade of the 19th. century, physicists had agreed that 
light was an electromagnetic wave. And out of the experi- 
ence with waves in other contexts had required a medium 
of propagation for it, the ether. Yet ether evaded all at- 
tempt of detection no matter how ingeniously the detec- 
tion devices were designed. On the other hand the notion 
of medium of propagation seemed so obvious to every- 
one, that some of the brightest minds of the time invoked 
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theories that would have kept the ether, but in a con- 
cealed form. The famous transformations of Lorentz were 
developed primarily to explain the null results of the ex- 
periments of Fiscau, Michclson, Michelson & Morely, etc., 
rather than as a foundation for the special relativity that 
emerged later. It was up to Einstein to think and later 
to speak out that our concepts have no legitimacy beyond 
what the complex of our experiences bestows upon them. If 
experiments do not reveal the ether it could be dispensed 
with. If experiments show that the speed of light is the 
same in all reference frames, that could be a fact to build 
the new physics upon it. 

In these early years of the 21st. Century, it would not 
be unwise to take Einstein's advice seriously. Accept the 
validity of newtonian dynamics in solar and similar sys- 
tems, where it has been tested, but look for alternatives 
in larger galactic and extragalactic scales, where it has 
not been verified. Many researchers have actually adopted 
such a point of view before. See Sanders and Mc Gauph 
(2002) for a review of Milgrom's MOND and earlier at- 
tempts. Here, we adopt a variational approach to the prob- 
lem. We amend the classical action by adding a term to 
it that eventually makes provision for the field energy to 
serve partially as its own source. The procedure is apt to 
a general relativistic generalization, in preparation for a 
future presentation. 

2. Variational formulation 

The system to be considered is a collection of point masses 
m,i at positions Xj, i = l,---,N; plus the field </>(x) at 
point x. Both newtonian equation of motion and equa- 
tion of gravitational field are derivable from the following 
action integral. 



J L[xi,0(x)] 



dt. 
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At point x the gravitational energy density is 
(87rG) _1 |V(/)(x)| 2 . By Einstein, this energy has an 
equivalent mass. We propose to entertain the possibility 
of using this mass or a fraction a of it as the source of 
extra gravitation on the particles. The effect on nij at 
position x l will be aGmi(87rGc 2 ) -1 1 V</>(x)| 2 |x — Xj| _1 . 
The contribution to the action integral from all space and 
all particles will be 



hnt = J L mt {xi,(f)(x.)} dt, 



L int [xi,<f>(x)] = 



|v<x> 



-d A x. 



(2) 



The dimensionless constant a is of the order of unity. It 
is introduced for possible later adjustments. Note that 
Gmi/c 2 is of the order of the Schwarschild radius of m^. 



Before proceeding further, let us emphasize that the 
argument given in composing 7j nt is for mnemonic pur- 
poses only. The formal statement of the assumption is the 
following: The classical action integral from which New- 
ton's laws are inferred is inadequate in galactic and ex- 
tragalactic scales. We propose to amend it by adding Ii nt 
to the classical action. The addition is a scalar, quadratic 
in field gradients, and contains the particle coordinates 
and masses. It will alter both the field equation and the 
equations of motion. 

The field equation for </>(x) is obtained by requiring 
the functional derivative of the total action, I + I int , with 
respect to (/>(x)to vanish. Thus, 
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By writing £(x — Xj) = 
integrates Eq. (3) into 
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Denoting the newtonian potential by 

irij 
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Eq.(4) integrates into 
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Equation of motion for rrij is obtained likewise, 
6(1 + I int ) 



6xi 



m,i 



_ , , . 1 a_ 
Xj - Vj0(xj) - 



|V0(x) 



-d A x 



= 0. (7) 



Newtonian limits in all cases are obtained by letting a = 0. 
To elucidate the significance of the proposed amendment 
to the classical action and the way it relates to mass dis- 
crepancy, several simple examples are worked out in sec- 
tion 3. 

3. Applications 

3.1. Gravitational field of spheres 

Consider a sphere of uniform density p, radius R, and total 
mass M: From Eqs. (5) and (6) the outside solution is 
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where s = aGM/c 2 is of the order of Schwarschild 
radius of M. For the inside, noting that <Pn{i~) = 

1 GM 



2 jp 



(3R 2 - r 2 ) , one obtains 



<f>(r) = — In 
dd> GM r 



dr 



R 2 R 



3 s 1 s r 2 
1 ~2R + 2RR 2 
3 s 1 s r 2 
~2R + 2RR 2 



r <R, 



r<R. 



(9) 



times larger than the classical ones. Having noted this, it 
is easy to write the first integrals of Eq. (11): 

Momentum : P = Mfx, = const. (12) 

Angular momentum : L = M? ff Xj xi; = const. (13) 
Energy : E = T + V = const., (14) 

T=\Y, M ^l V = -\GY J MfMf/n 



3.2. Motion of test satellites 



Virial theorem : 2T + V = 0, 



(15) 



Equation of motion is obtained by substituting Eqs. (8) 
& (9) in Eq. (7) and reducing it. The result is 



GmM 



mr 



r — s 



(10) 



where u(s/R) is derived, numerically calculated, and plot- 
ted in the Appendix. As the central body contracts and 
R reduces from an initially large value to its allowed min- 
imum ^s;u(s/R) increases from zero to 4. The first term ^ Concluding remarks 



where 'bars' on T and V indicate averages over system 
members or over time periods. 

Increased effective mass in Eq. (11) and the virial the- 
orem may prove useful in explaining the missing mass or 
missing gravity issues, in the velocity curves of galaxies 
and/or velocity dispersions in clusters of galaxies. 



on the right hand side of Eq. (10) is recognized as the clas- 
sical expression. The second term, u(s/R), is the pivotal 
one and as mentioned it could be four times as large as 
the classical term. In newtonian parlance this amounts to 
saying that a compact object may, gravitationally, present 
itself up to five times more massive than the same object 
in a diffuse state. This is what we propose as a partial 
answer to mass discrepancy! 

The third term in Eq. (10) arises from deviations of the 
force law from r~ 2 . It is infinitesimal compared with the 
previous two terms, though it may observably perturb el- 
liptical orbits in secular time scales. 

3.3. N-body problem 

Consider a collection of spheres each of mass Mi, radius 
Ri, and positions Xj. Assume the system members are well 



separated such that |x$ — x. 



> Ri for all i ^ j. It 



is demonstrated in the Appendix, that on account of this 
assumption, the integral appearing in Eq. (7) splits into a 
sum given by Eq. (A. 7). Thus, one finds 



x i -G^Mff^4 

^-^ J Xj — y|j 



X, 



0. 



(11) 



where Mf = Mj(l + uj),Uj = u(sj/Rj). Actually a term 
of the order G J2j Mj Sj | — s } ; | ~ 1 , is neglected in Eq ( 1 1 ) , 
a) because of the extra small factor Sj and b) because of 
its steeper, r -3 , decrease with distance compared with r~ 2 
dependence of the retained terms. 

Albeit the approximations, Eq. (11) is identical with 
the newtonian many body equation of motion, except 
that instead of the classical masses, the effective masses, 
Mj S = Mj(l + Uj), play the role. Depending on the de- 
gree of compactness of the constituent bodies, the effective 
masses could, in the case of uniform densities, be 1 to 5 



... it is contrary to the mode of thinking in science to 
conceive of some thing that can act itself, but which 
cannot be acted upon. Albert Einstein, 1922. 

The quotation is from Einstein's argument to set the 
stage for the rejection of absolute space-time continuum 
of special relativity in the presence of matter, and his 
admiration of Mach's reservation to ascribe an absolute 
meaning to inertia irrespective of matter elsewhere. In an 
attempt to amend newtonian dynamics, we have followed 
Einstein's advice and tried to alter the laws of gravity 
and of motion simultaneously. The interaction lagrangian 
of Eq. (2) contains both the particle coordinates, Xj, and 
the field gradient, V<j>. 

Noteworthy in L int is the field energy that eventually 
appears as the increased effective mass, M cS = M(l + u 
in Eqs. (10-15). We recall that the main issues in the 
analysis of the velocity curves of spirals are: a) insuffi- 
ciency of the observed and/or estimated stellar masses 
to provide the required dynamical effect, and b) the flat- 
ness of the velocity curves beyond what the distribution 
of observed masses would permit. Similarly, in clusters of 
galaxies, again the observed masses are much smaller than 
the virial ones, inferred from the velocity dispersion in the 
cluster. The increased effective mass discussed above could 
provide a partial answer in both cases of luminous mass 
deficiency. The flatness of the rotation curves is a mat- 
ter of the distribution of effective masses. For example, if 
there are enough low luminosity collapsed faint objects in 
the outskirts of spiral arms the velocity curves could be 
flattened enough. 

The field in Eq. (6), unlike the Newtonian one, is not 
proportional to the mass of the constituent members and 
does not admit of superposition principle. Even in the 
case of approximate superposition in systems with widely 
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separated components, the field of each component is a 
highly nonlinear function of its mass. This nonlinearity 
could prove useful in discussions of Tully-Fisher relation 
in spirals, an empirical relation between the brightness of 
the galaxy and the asymptotic orbital speed at the flat 
end of the velocity curve. The Tully-Fisher relation, with 
the assumption of a further empirical mass to light ra- 
tio for the galaxy, translates into a power law between 



the asymptotic speed and the total mass, v a 



oc 



2.5 < (3 < 3.5. In the present dynamics, (3 will sensitiv- 
ity depend on the compaction factor, the fraction, and 
the distribution of the compact members of the galaxy. 
Same considerations and reservations holds for the Faber- 
Jackson relation, another empirical relation between the 
velocity dispersion in cluster of galaxies and the luminos- 
ity of the cluster. 

In section 3.3, conservation laws were derived for the 
case of approximate superposition. One can, however, do 
better. The total lagrangian, L + Li nt , is invariant under 
time translations, coordinate translations, and coordinate 
rotations. Therefore, seven first integrals of motion, cor- 
responding to energy, momenta, and angular momenta, 
should exist. This point of view will be presented else- 
where. 

Birkhof's theorem states that in newtonian and rela- 
tivistic regimes the gravitational field of a spherically sym- 
metric system is independent from the internal structure 
of the system. This is not the case in the present theory 
as it is highlighted by the presence of u(s/R) in Eq. (10). 
This may open the possibility of studying the history of 
contraction or collapse of an object by logging the external 
gravitational field of the object as it evolves. 

From Eq.(7), the effective force on a test particle is 
derived from a potential. It could be eliminated in a freely 
falling frame of reference, meaning that the dynamics 
developed here admits of equivalence principle. This is 
hardly surprising. For, in the action integrals of Eqs. (1) 
and (2) one single mass m, is used to compose the kinetic 
and the gravitational potential energies of the system. 

Finally a criticism: a) The mass equivalent of energy 
is a relativistic concept. Yet, a gravitational effect was at- 
tributed to it through the law of gravity of Newton, b) The 
effective mass attains its full significance when the source 
object shrinks into sizes of the order of Schwarzschild ra- 
dius. Logically, a covariant general relativistic approach 
should be adopted. Presently, we arc looking for a field 
equation and a geodesic one in which some sort of mu- 
tual interactions between distant points, similar to one in 
Li n t are taken into account. As in Eq. (2), there will be a 
free parameter a in the formulation. We will require the 
emerging dynamics to reduce a) to the general relativistic 
one in the limit a — + 0; b) to the present dynamics in the 
limit of weak fields but a/0, and c) to the newtonian 
one in the limit of zero a and weak field. This will be 
presented elsewhere. 



Appendix A: Reduction of Li nt , Equation(2) 

We have already seen the formal role of the interaction 
lagrangian, L int in altering the equations of motion and 
of the field. Its numerical values are needed in the study 
of orbits, Eq.(7); virial theorem, Eq. (15), etc. Here, we 
calculate it a) for a single uniform spherical mass, and b) 
for a collection of unform spheres. 

a) For a spherically symmetric field the non newtonian 
integral appearing in Eq. (7) is 



; n - — — 
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-d 3 r'. 
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In an expansion of [r — r'] -1 in Legendre polynomials, 
only the I — term will have non vanishing contribution; 
for d<j)/dr' has no directional dependence. We carry out 
angular integrations and split the radial range of the in- 
tegration into three intervals < r' < R, R < r' < r and 
r < r' < oo. Thus, 
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The expression for the field gradient in the first integral 
is the interior solution of Eq. (9) and in the other two are 
the outer solution of Eq. (8). The remaining mathemat- 
ical manipulations are elementary, though lengthy. One 
eventually arrives at 



GM 
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where 



u{y) = -- 



1-4/— — — arctan ' 
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(A.3) 



(A.4) 



The parameter s/R indicates the degree of compaction 
of the central body. As the body contracts from an 
initially dispersed state, (s/R) ~ 0, to its maximum 
allowable compact state, (s/R = 2/3), u(s/R) increases 
from zero to 4. A plot of u(s/R) is given in figure below. 
In newtonian language, the effect on the motion of an 
orbiting satellite is to make the central body look (1 + u) 
times more massive than what it docs in newtonian 
regime. 

b)Multi-component systems: In Eq. (A.l) The main 
contribution to the integral comes from the immediate 
vicinity of the central body. As one approaches the grav- 
itating body from outside, \d(j)(r') / dr'i 2 grows roughly as 
r'~ 4 up to the surface of the body and fades out to zero 
at its center. Assume a system consisting of two well sep- 
arated spheres (Mi,Ri) at X! and (M 2 ,i?2) at x 2 with 
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Fig. A.l. u(s/R) versus s/R, the degree of compactness 
of the central body. 

|xi — X2I 3> i?i&i?2- From Eqs. (5) and (6), we read 



|V0(x')| 2 




Again, the main contributions come from the vicinities 
of the two spheres. In the immediate vicinity of sphere 
1, however, G 2 M 2 \x' — xi|~ 4 is large. While, with 
the assumption of large separation of the two spheres, 
G 2 M| |x' — X2 1 ~ 4 is insignificant. Vice versa for the vicinity 
of sphere 2. The expression containing (x' — xi).(x' — X2) 
practically adds up to zero upon integration over x', 
because of its directional dependencies. We conclude that 



GM x ui GM 2 u 2 ( Si\ 

knt{x)~- - + - r, Ui = U[—\. (A.6) 

|x — xi| |x-x 2 | \Ri J 

Similarly, for a system of many, but well separated, com- 
ponents one obtains 

M*)«g£ u *(*W (a. 7) 
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